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AN AXISYMMETRIC, NUMERICAL MODEL FOR 
A NON-HYDROSTATIC BOUSSINESQ OCEAN 


1. GOVERNING EQUATIONS 
The governing equations of the axisymmetric, non- 
hydrostatic, Boussinesq ocean model are 
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where VY" = Re Siena other symbols are listed in Appendix \. 
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Above, the density anomaly b is defined according to 
see ee 
r 
b = —————- g , pls) 
Po 
where o,(2) is a reference density and is a function of 
depth only. Brunt-Vadisilla frequency N- is defined as 
(1-6) 
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The continuity equation is that of the incompressible fluid, 
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2. THE MODEL GRID 

It is determined that a tully stagevered grid is most 
expedient for storage economy fora yviven spatial resolution. 
As shown in Fig. 1, the radial (u) and the tangential (v) 
velocities are defined at cross points, vertical velocities 
(w) are defined at open circle points, and the pressures (p) 
and density anomalies (b) are defined at blackened dot points. 
This grid system has the following advantages: 


a) it saves storage for agiven spatial resolution 


Se 


b) it is very economical in terms of number of 
computational operations for the finite 
ditference (FD) equations of (1-1} to (1-4). 

c) it is very easy to specify the boundary 
conditions, 

d) the pressure diagnostic equation, of the 
elliptic type, can be reduced to the standard 
rorm, and 

©) there is no spatial separation of solutions 
on the grid. 

In order to consistently index the grid points, we 
let index pair (ij) represent the i-th point in the 
r-direction and j-th point in the z-direction. In addition, 
m is the maximum number of points intner-direction, and n, 
the maximum number of points inthe cz-direction. Therefore 
there are m x (n-l) points for radial and tangential 


velocities, (m-1) x n points for vertical velocities, and 


(m-1) x (n-1) points tor mass distribution (b and p). 
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3. THE FINITE DIFFERENCE EQUATIONS 


The leaptrog, or centered-in-time, integration 


scheme for the inviscid terms and the forward-in-time 


integration scheme for the viscous terms are used. The 


scheme is described as 
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A second order, or centered-in-space, 


scheme is applied to | 


derive the tendencies in (3-1) according to (1-1) ~ (1-4). 


(a) The Equation or Motion in r-direction 
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(b) The Equation of Motion in 9-direction 
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where 


The Equation of Motion in z-direction 
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4. DERIVATION OF THE DIAGNOSTIC EQUATION FOR PRESSURE 

The nonhydrostatic pressure at time t is needed to 
compute the pressure gradient forces in (5-2) and (3-5). 
To "recover" the pressure from the motion fields, we make 
use of the continuity equation by differentiating (1-7) 
with time we get 
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which can be written in rinite difference torm for a mass 


point ij as 
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Substituting (5-2), (3-5) and (4-2) into (4-1), we have 
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"we obtain the standard fgorm of an elliptic equation in FD orm 
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Equation (4-5) can be solved numerically by the SEVP solver 
iMadala, 1973), providing the boundary conditions are 
properly posed. 

The conditions for the four boundaries are determined 


according to the following assumptions: 


(a) At (Ti); = (ray = 0, the natural condition for the 
cylindrical coordinates calls for u = v = 0 = 3u/3t = 3v/9dt, 
the gradient balance requires that (oP/ar) oy = 0. Therefore 


an extra column of P is needed 
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ib) At helene (Today assuming both the horizontal 
divergence and the vorticity are continuous, i.e., + L a = 0) 
aX oLlosvr : 

and = = a = 0. These lead to 
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Note that if by is set equal to cero, (4-7) describes a 


non-divergent and csero-vorticity boundary condition at 


r= (Py dine Once Vii is determined, a gradient balance at 


c Coe requires 
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where a column of dummy points Pini has been introduced for 
computational purposes. The second part of the RHS of (4-8) 


is thus the forcing function at (Tye)in tor the elliptic 


equation (4-5). 
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Following the same deduction between (4-1) and (4-5), we get 
an expression similar to (4-5) with the second team on the 
LHS and Oey in the RHS absent. Thus, Bag can be obtained 
by the same SEVP solver by setting Cr) = 0 and Oia = 0. 
(d) At top We za/st Wi * 0. Following the same line 


of reasoning as in (c¢), we obtain ae by solving (4-51 with 
a Pes = 0) and G.. = 0. 
e n intl Yin 

In summary, the elliptic pressure diagnostic equa- 


tion (4-4) is to be solved with the following boundary 


conditions 


1) \t r = 0 aA = ra 1.e., (4-6) 
2) At r = Lai Pini = Pp oe + function (Vin) {4-8) 
3) At 2 = 0 AZ, = 0 and Gay = 0 
4) At 2 = (2), Ce = 0 and or = () 
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LIST OF SYMBOLS 


an array of constants, varving only in r-direction, 

defined by (4-4), used in (4-5) 

an array of constants, varying only in z-direction 

defined by (4-4), used in (4-5) 

an array of constants related to ry and AY, used in 

(4-2) 

an array of constants, used in (4-5) 

density anomalies, defined in (1-5), cms 
\ 

an array of constants, varving only in r-direction, 

defined by (4-4), used in (4-5) 

an array of constants, varying only in z-direction, 

defined bv (4-4), used in (4-5) 

an array of constants, related to ry and Ary» used 

in (4-5) 

Coriolis parameter, st 

gravitational acceleration, cm 57° 

an index, denoting i-th point in r-direction 

an index, denoting j-th point in z-direction 


> - 
horizontal diffusion coefficient, cm™ s : 


> os 
vertical diffusion coefficient, cm” s 
left hand side 


the maximum number of grid points in r-direction, 


upper bound of i 


15 


1 

m, m-2 

Ne Brunt-Vadisilld frequency, s 

n the maximum number of yvrid points in z-direction, 
upper bound of j 

ny n-l 

n, n-2 

7 a? 

p pressure, dvne cm ~ 

RHS right hand side 

Tr radius, cm 

rs radii of momentum points, cm 

T> radii of mass points, cm 

Ary distance between two horizontally adjacent momentum 
points, cm 

Ar, distance between two horizontally adjacent mass 
points, cm 

SEVP stabilized error vector propagation 

t time, s 

At time interval, s 

u radial velocity, cm Cie 

Vv tangential velocity, cm 3! 

Ww vertical velocity, em a 

z height from ocean bottom, cm 

zy heights of circle points, cm 

Z4 heights of cross and dot points, cm 
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haar 


distance between two vertically adjacent circle 


points, cm 


distance between two vertically adjacent cross or 


dot points 
density, gy cm 
a constant density, 1 g cm 


a reference density, varving 
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APPENDIX A — FORTRAN CODE FOR THE NON—HYDROSTATIC MODEL 


A listing of FORTRAN code of the ocean model. The 


major functions of the main program and subroutines are as 


follows: 


OCEAN 
INIT 


START 
| PUTOUT 
} MAP 


ADVECT 


DIFF 
PRESS 


MATINV 
BSM1 

BSM 
‘BSM 


4 


to 


ua 


FRWRD 
BOUNDV 


CHECK 


main program, calls all subroutines, manages 
job flow, controls input/output. 


sets up independent variables, defines con- 
stants 


defines initial conditions 
gets various fields ready for output 
prints 


computes all inviscous terms, except for the 


pressure gradient forces 
computes horizontal and vertical diffusions 


solves the pressure diagnostic equations and 
computes the pressure gradient forces, appears 
only in the non-hvdrostatic version 


inverts matrices 


Used in SEVP method 


matches forward 
sets outer boundary conditions for momentum 


checks if the time step is linearly stable. 
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APPENDIX B — FORTRAN CODE FOR THE HYDROSTATIC MODEL 


The hydrostatic version of the model can be obtained 
by simplifving the non-hydrostatic version. In the hvdro- 
Static version, the equation of motion in z-direction ;1-5) 
is reduced to the hydrostatic equation 


-i@-=-p .B-1) 
eee 


Instead of solving the elliptic equation (4-5), the pressure 
p is thus obtainable by vertically integrating (B-1}. The 
vertical velocity w can also be computed by vertically inte- 
yrating the continuity equation (1-7). 


The FD forms of (B-1) and (1-7) are, respectively, 


t 
Pehl 


Ty 
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CALL INIT 
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MORIZANTOL ACVECTION FOR RADIAL VELSCITY 


O* to Jat,nt 
O53 lq YeQ, 
LO VRYCT eI) Beg 25a c (VRC lL, JI oVACToged) DeCVRAC Tp JIOVA (Lote JS) ,CAL(Tot) 
1 OOVR Tate JIOVACT JI IMCVRC Toye JIOVA CY, JII“ZORI(T)) 
2 oVR3(T J) 


WQRITZINTAL AOVECTIIN FOR TANGENTIAL VELOCITY 


ce 29 Jat,Ni 
VP 20 JTs2,mt 


20 VIMT eT Bede2Sac (VR (Te SI OVACTess J OCVT CT ed IoVT(Leted)) sOP I Clot) 
1 SCVACT Ons JISVACT JIIMCVICT ote SIOVTCT JI I/O%N (LD) 
2 ovTscr,d) 


MARTZONTAL AMVECTIAN EAR VERTICAL VELSCITY 


O° 30 Jae,nt 
D8 30 382,42 
SO VZ9CT oJ Meg e2Se( (VRC eo JI ovACIe Jory eCVZCLedov2Cleted)) sCRACT) 
1 OCVACT ON SIOVACT OL etapa cvZCTetetiov2 (led) /0K2S¢0I 41>) 
2 *vZ23(1,J) 
DD 40 gaz nt 
49 OPE RM SE Uee eet ety Teas eeed ney imtnn brane we iendaecee ta? 
ov2 31d) 
ae Se aden 
Sa eee Jy Bad 25ae(VR(MT J ouF (Mee Jory ecvZCMl Soul (M2,J)) /OR2("1) 
ovedenn, J) 


HORIZONTAL ACVECTIOAN FAR BUBYANCY 


of tu Ja1,%t 
CO 60 Isa,N2 
60 BSCIASIBASCT SieO FACVACT  SIeCBC Le Joh (Tet, J) NRE CT) 
! PVR (Tope J MCR CT Or, JISACT, JI I/0R2CI41)) 
Of Tu JaieNi 
70 BCMA IMATCI, SI oO SAVE (2S, JIACRC2, JBC ied) IZOR|A(2) 
ce 40 Jai,ni 
Ry BCMA DIBASCML JOO eSPVR CMI gS (OCH J I@BCM2,3) sCR2(H1) 


VERTICAL ACVECTIAN FAR RACTAL VELSCITY 


CA 90 Jayne 
ne ay Ise,™1 
9g VRSCT oD OVAR CTs edeSSal(VA(TelsJIoVICTeJIIMCVRCT > JI@VArE,Jeoi)) 
! AOZACIVOCVACT, Jol yovZ(Laot Jol) ecvACl Jot yeva(y,s)) 
2 MNZ2C JON)? 
Qe 9g fag,at 
5 VRECTpPIBVAS (CT 1d eSaCvZCLe2revzZ (let, 2) a(vACl, 2 ,evA(t41)) 
7022¢2) 
ar Fo [ag,nt 
Ge VASCT ANI BYR ECT NSEC eS aCVZC TOL NI JOVZCL oN) DOCVACT SL yovACEeN2)) 
1 smzZ2acne) 


VERTICAL ACVECTION FAR TANGENTIAL VELOCITY 


31 


0001000 
0002000 
0003000 
0004000 
0005000 
0000000 
0007000 
0008000 
oovecag 
0010000 
0031000 
0012000 
0033000 
0034000 
0035000 
0030000 
0027000 
0618000 
9919000 
0020000 
0021000 
0o22000 
0023000 
0026000 
0025000 
0020000 
0g27000 
0928000 
0029000 
0930000 
oodi1000 
ood2cce 
0033000 
oosacee 
0055006 
0036000 
0057000 
0058000 
0959000 
0040000 
94 {000 
ogezcge 
Oo4yood 
oa4a0ce 
09045000 
0046000 
0047000 
0046000 
0099000 
0990000 
0091000 
oos200c 
0093000 
99940090 
00999000 
0096000 
0097000 
0036000 
0039000 
0000000 
0001000 
962090 
0983000 
99ea0G0 
00050G0 
0006000 
0007000 
0008000 
0069000 
0070000 
0071800 


CARO [mack FILE ECTTOR(CTFEH) es VEISTON 05,29 OHATERIOs26/42 


eae “EMBER ACVECT 


Tae 
730 
748 
7$a 
Joe 
TTa 
780 
790 
AQe 
ele 
@2e 
Aye 
aue 
Se 
aoe 
a7" 
eRe 
ede 
Q0a 
Qie ¢ 
92e C¢ 
930 C 
que 
95a 
Gee 
G7e 
abe 
99a 
1n0e 
lgise 
Inde ¢ 
tO3e ¢C 
194e ¢ 
to5e 
Joon 
love 
1940 
109e 
1100 
Vsie 
li2e 
L13e 
1140 
1156 
t1ee 
Lite 
118e 
L192 
tade 
tate 
1ade 


ann 


RAN 


aac 


sae “EMMEN AQVECT 


Le tyd 
£8 1y9 


Ja2,2 
Ts2,%) 


TIPEalasissi gst? 


LOG VT ECL oe TIBVTSC  JeOeASal (Villers JovZ(LsJIICVTCT oe JIevT(1,Jel)) 
1 ANZTSNCIIOCVICT JeotvevE (let, Jor) relvrilsdotpevietge dd} 


é 
g* 105 


AOZYCJot)) 
1a2,%1 


CUS VISCL ot davrs( I tpeGedSalvZCle2ovZ (lel, 2) acve(L,2yeowrit,1)) 


1 
C* lige 


PVG VISCTONDBVTSCT Nd Der e25aCVZCT yh OVZC TOL eNS) IOC VICI NE pevr(l,N2)) 


mpz2(2) 
Ts2,"1 


enZgeny) 


VERTICAL ADVECTION FQR VERTICAL VELACITY 


C7 tau Jeet 


om tad 


Tat,“4 


V4H VZZCTODIRVZS (CT Je VeSSeCCVIECT Jer vovecie SI eculcLesievZ cL Jot)? 
JOEY CImHVOCVACT, Jos vov2Z(T ,Jrecvicl Jet yevi(t,J)) 


1 
2 


ca 394 
cs ise 


enZ4039) 


VERTICAL S&OVECTION FAP 8 


tea, 2 


Tat,“4 


WSQ AVCT AD ISASCT, Fea SeCvZC Ls JreCBCLeyo@ (lL, Jelyys0Z2¢s) 


i 
Te Lyn 


OV ZT Jot yaCAl( pr Joi oA(l JJ I/CZ2CSe3)) 


Tai."4 


Ted AICLo»VBOB(T Lou, SevZ(L,2a(Rallsgieracts rip /022¢0g) 


Ecco i AT] 


rst,“ 


VTE CBC AIBA CT MA eMeSevECT NL IMCACT Mi enlLeN2) I/O Z2(h1) 


INERTIA TERMS FER KORIZGNTAL MOMENTUM 


GM 110 Jatevt 


ce 140 


VASCEe JIMVRICT JOVI CLS aCVTC LoS sRi CL) oCORT) 
tto VIPCL eS avr S(T, JeVACT, J aCVICT se SIAL CED OCORI) 


Ta2,“1 


ALOVENCY TERM FAR VERTICAL ACCELERATION 


cf 120 Jad,*1 


of teu 


Tat,*t 


V20 V2FCTo SI avZ3( Ld Oe SACP CL J oBlIeJeor)) 


oe 130 
of au 


SThatyeEtCavy4n TERM 


Jei,'t 
Tat,“1 


13d ATCT SIMBSCI J ea SMCVZCT eI) OBVES CTP OVA(T Sos Oval So1)) 


RETURN 
ENO 


a0DEQ TA 


seuece 


122 AECORDS 
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0072090 
0073090 
oarddgd 
0079990 
CO7Tecoe 
aorroge 
0078090 
0079090 
oaeonad 
0081000 
cor2000 
0083000 
0064000 
agegagd 
ogeecod 
0067090 
ogeeced 
90899900 
0940090 
0093000 
oov2og0 
covsuad 
ogwergd 
09960900 
0046000 
0n47Cg0 
0ovecod 
0099000 
0100000 
v19f0gd 
arc|ecge 
g1U3000 
O1G40ge 
0109000 
a10@c00 
0107900 
0104000 
0109000 
0130000 
0111000 
Or1sgcge 
0133000 
91146090 
o121§c00 
0136000 
9117090 
0136000 
0119090 
oré2geae 
0141000 
0142000 


CARD LmaGle FILE EXTTAR(CLFEN) eo VE2S]AR 05,29 OATERIOs20/02 TimEatagigsi289° 


aoe MEMBEX PUTOLT 


1s SL ARGUTINEG PLTSLT 0001090 
2a PARAMETEN MB2i,n a6! ogv2000 
3a PARAMETER Mistoy ,h2EMe2,NISNol ,N2SNeg 0003000 
ue COMMANSINESVRI CHAT IOVILCM NES VZUCML aN) oBLCMESNEDVRAC Hot I )y 0004000 
Se 1 VT SCMND eVZ2CMI AY) eBSCML ONT GVASCMANL) oVTS(M NY Dy oouscod 
oe 2 VZZCME AR OIECMIGNED SPOMIANI) 6000090 
Te COMMON THA ZAE(M) h2C41) DRE CME) DRAM) ZICN) pZ2CN1) OZICN1) gOZ2CN) 0007000 
Ba COMMONS THR ZRH BRAREND ) pAV2CN) ALPHA, ENDAGONDR CORT sGyhk(™),2K(N) 00U8000 
Cry CAMMONJZENP JOEL T ATIMESITIMES ISTEP, TSMO, I TAPE, TEV 0999000 
{Oe €¢ 0010000 
tte ¢ TeT$ SUBROUTINE PRINT GLT FIELDS FOR A GUICK LOBK 00110900 
tea c 0912000 
13e CIMENSIAON INL (A) 0013000 
140 IAG EMRMATCAS/4—¢" FACIAL VELICITY (OMsS) aT Ta',loe' wt) 0014000 
15 79S FIGMATC////," TARGESTI OL VELMCITY (C*%sS) AT Tet,los! bt) 0015000 
1oe Pig FARMATCAS//5" VERTICAL VEL“CITY (CMsS) AT Tw’, ley! wt) 0016000 
17a 715 FVOMATC///7_et ELEYANCY FIELD (9,001) at Te, Toe! wt) 0037000 
168 725 ENRMATCS/S/ 9" ERE SSLRE (oso DYMEsCHneg) AT TIMES’, 1o,!' wt) 00tacge 
149 720 FVAMATCIFIG/SS/, aC he todesorete ALTOLT AT TIME Biele,! ge, 0019000 
ece 1 Fagdet Cav ISte® a',I7,! osereseoooes!) 00¢€0000 
ete ULAVaXTIMEsHOAND LoL OCC1 Ood100C 
220 PRINT T20,TTIME, Cav, ISTEP 00e2000 
23a CS to Jei,vt 0923000 
ean 09 Ju Tet, Ogeacyd 
2Se 19 1DUMCl,J avI2(T,J) 0o0e$00c 
260 PRINT TOO,ITINE O00¢etod 
27a CALL PAPCINtImsay, 22> oN) 0027006 
26e O7 2u Jet,rt 0gestoo 
29a no 2u late 9029900 
JOe 2g POuM(l-Jdavraci.J) 0050000 
gis PETtT PISgITIVe 0031000 
320 CALL MAPCTOUM SE Eyl ag yh) 0052000 
330 L° 39 Sete 9053000 
W4e m* Sy Lebent ooSucge | 
350 Bo TMM CL, pav7Qtt, J) 0035%00 
toe OF TUT Mus ll EME 0036000 
37a CALL MAP CIN wm, CZ elie“ 1%) ogs7co0 
yAa ch 4 SEL gt °998000 
398 no dy [aie ¢n39000 
afte ag (CuCl J shsCTJieret? vosucoe 
Pa) PRINT TIS, IT YE oosrcno 
uae CALL MAPCIOU Ae eZ2eM1 Nt) oo#ecoo 
age mY Sy JRr— il 60943000 
ada o* 59 [er,st Coeucge 
use Seq usM(l, Josef re Jye1,Ee! 9945000 
vor PRINT T25,I1TIME 9046000 
ave CALL MAP CINLMyRe,2ae “tpt? nga70g0 
ube RETURN 104acod 
ada ENC oosetag 
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CARD LMAGE FILE EF TTAR(CIFEK) es VESSTAN 95,29 


eae MEMEER “ae 


ta 
de 
Sea 
Ge 
Se 
6a 
Pe 
Be 
Ga 
10" 
1le 
tas 
130 
ida 
San 
Loe 
t7e 
1608 
190 


3o 


GUAROLTINE “a Cag gle ™4 NN} 
PARAMETEN MBQ1,K84) 
OQIMENSTAN RCM), 20NN) 
INTEGER & (M,N), TECMI,TZ(N) 
FARMATCImMS,7xe2815) 

FARMAN (inS,1u,3x,2515) 

MP g*}hO(25,%™) 

Of to Tay,nP 

TRL IMAC LIel Eager! 

Of 20 Jei,Nr 
TZ(JISACII #1 Ewe o,t 

PRINT Fu 

PRINT TO,CIPCT),1ege"F) 
PRINT 7C 

Of 30 JJat,*n 

JOnNeleJdJ 

PIUIMT BOLIZC SI, CAC Tod), Tai,*P) 
RETURN 

Lael) 


34 


O*TEwtusges62 


TIMEStagsgri2se? 


ogo1led 
09002090 
8003000 
0004000 
0gugcad 
0006000 
gquregd 
1008000 
0009000 
eo!00g0 
Co11dg0 
2012000 
Og1 3000 
ooluvot 
0015090 
0016000 
9017006 
0918000 
0019090 


CARO ITHaGE FILE ECTTAR(CTFEM) em VEPSTAR CS.29 HATESLUszesba TIMEStusizer2s9? 


tee PENMEN PRESS 


1a 

ae 

3a 

ae 

Sa 

6a 

Ve 

Re 

Ga 
100 
ite 
120 
130 
14a 
1S¢ 
168 
ive 
18e 
198 
20e 
aie 
eee 
ese 
24e 
2Se 
~e6e 
ave 
28e 
aoe 
300 
Sie 
320 
330 
yae 
35a 
Joa 
J? 
yaa 
390 
ada 
Gia 
uae 
a3a 
ude 
aSa 
abe 
ata 
ufe 
49a 
SQe 
Ste 


aaAna 


Nanna 


Cc 
c 
C 
C 


aan 


SEBROL TINE PRESS oovrced 
ecu2rod 

TeTS SUSRALTINE SETS LP FRWCING FUACTISAS bA0 FOUNG AR YOQUSC UG 

CAADITIONS FAR ThE PRESSURE DIAGNASTIC EGuaTLaNsS FAR Sd0t4C00 


005090 
PORAMLYER MB2t,K BZ} COvet00 
PARAMETER Miao, waered,Niatet,23nez o0u7C00 
PAHAME TER MAL K aS, AULA I SNhL KO] oguaeccd 
Reabet RCOW,RINVRINVI,HTILCA,UL4Y1 Ceusgtgo 
CAMMONFZELESVRE(H AED OVTTCMAS EDS VZICMP ON) SALOME NT), VR2( My H), Ouroccd 
1 VTQACM NT eVZZCMI ND SBSCML ONT Do VRS(MINI) gp VTS(M NID, oor100d 
2 VZGCML Neb ECE ONL) POMS DE) 0012000 


COMMINAT HISAR (My RAC%1) pCOLCML)D ,DOPSCMISZLCN) ,ZZ2CN1 oA CZICN1),0Z2°0N) 2013000 
CUMMONSTHH AER  hRER (NT) AVACN) ALPMASBNOASGNOR CARE GRE CM) ,ZKCN) 0010006 


0015600 

»PenP [S$ THE SIZE 8F X AND F 00346000 

PPEMLOE, “PEN ouhL Kol 0087000 

0018000 

PARAMETER “Pg%o1,AOBne} 0039000 
PARAMETED “Oy ant of, "Pcm*Pes, KP LENPal,'P2gNPad 00¢0tQ0 
CRMMANZEVOSQT ay (hE SMP QQNALK) pRAINVI CMP 2 MP2,NHLK1) »PCSR (MPSS), 0o¢1000 

1 STILE A CMOS), F (MP SNP) NASTZ2(NBLK),IS(CNELX) pSLHF(NBLK), 006¢2000 

2 TE CNELE ) FIL CMP) pFINCMP) F221 (NP) SFAM(NP) SAK (MR) ,AY(NP), 1023000 

3 PACHE AB) CXCME) CV (NP) 0oe4u0ce 
OIMENSTIAN DUMMY) (PPA ,YP2), (MPN) 00¢5900 
ETLIVALEACE CoUMM YT, RIAV( 1,1 ¢NALKD) 026000 
PARAMETER HPMP gy DamPp 0ue7000 
CATA DUMMY IL /MPME AG spt sMPMP OU LY 0028000 
Cara NCALL 07 0029000 
"ASTZ2 REPFESE"TS KAUMBER RF INTERTAR GRIC POINTS IN EACH BLOCK IN XeDI0050000 
N2 REPRESENTS MUMRER CF INTERTOP GRID POINTS IN YeDIRECTIAN oosicee 
NALK FEPRESENTS NUMBER OF BLACKS IN XeOIRECTIAN 0032000 
ThE VARIABLES Alt, Ath, A21,a2M TAKES THE VALUE 9 FAR OIRICKLET 8,C, 0033000 
NO 1 FQN NEUMANN &£4C, Sf THEIR RESPECTIVE BOUNDARIES aty CORRESPOND00S4000 
Jat 41% TA Jan a2t tT) Tal azn =o r@ Tee 0035000 
BOUNDARY CONCITIONS ARE 0056000 
NCP op dmCheatpsax( r,s oA ge(eCl, 2 ,eF 19 (2)) 0037000 
ACLotP BC LeALY aX (To hP Atha (X(T NPOt)OFIN(T)) 0038000 
RCVo J aCteomar ran (r sass a(x (2, Jor 25 (5)) 0959000 
X(MP SBC 1eAQY) OK (MP, J) OAQMACK(MPaL, JI OFAM(J)) ooscced 
NCALLANCALL 41 0041000 
0g42ng0 

DEFINE THE FRPCING FUNCTION SF THE ELLIPTTC EGUATIBN 0043000 

0944000 

NF 10 Jat,nhe 094§000 

cf to Te1,482 0046000 

t¢ Pik Ala aia ie ol gs rh be an er ane ee Jerys021¢9) 0047000 
PAXCropeGRaCLaVR3CL JpevZ3(LeJI/0Z1 (J) aRne 6048000 

0049000 

SET LP AN INITIAL GUESS 0950000 

09051000 


35 


4 


eee MEMHER PRESS 


Ses 
S3e 
Sue 
$$e 
Soa 
Sve 
§ae 
698 
ede 
eis 
eze 
e3a 
e4e 
o5a 
66 
e7s 
o8n 
ede 
708 
vie 
ree 
730 
740 
750 
740 
174 
7688 
780 
e0e 
ata 
@2e 
e3e 
aaa 
aSe 
46a 
a7a 
abe 
eon 
908 
gta 
G2e 
93a 
Cue 
S$e 
968 
G?e 
qAe 
a9" 
1c0e 
icte 
luds 
{ois 
16d" 
icSe 
1cos 
1078 


anan 


aan 


CARU IMAGE PILE ET LTTOR(CIFE™) we VERSTAA 05,26 


20 
Jo 


005 


191 


i106 
“4e 


t3e 


CATERI0/26782 


TFCNCALL GT,tIGR Te 3¢ 
C3 20 JBL nee 

nA 29 yar,4Fa 

CTL eJeorssP (Ig) 
CONTINUE 

CALL ZILCHCF LL RF) 
Care ZILCH (EIN, &P) 
Cath ZILCHCF21,KP) 


OEFINE THE FORCING AY BCUNDARY SO TeAT TreEkE 
IS GRACIENT BALANCE aT OUTER SOUNCARY 


C7 605 Je2,"P1 

Fem C Jahn MaNRa ch pavl]acMsJetyaC(vracmytetysRi(m>oCORy) 
CANTINLE 

4118! 

AINSI 

a2isi 

aAgrsal 

27 VOL Sm aenOt 

BP C2,d ISAACS, J) oan(2)eael 

F(QeJI OF (2, J poanceeFet( sada 
MCteJIB(Loeaat yan (ted) 

ABQ(*Pey, J BARNET, JOCK (“Pel )edam 
ECMPol, J) BF (MPa, JSeCU( “Pet aFaX(sjyeage 
£CMP SSCL Oe om) an (MP, J) 

CINTINVE 

°° 102 T82, “Ph 

PRCT, 2)BHACT, 2) e4v (2) eat] 

FC Se P2IaF CT, 2p ear (ereFiir(rjyeatt 
RCLotdaCisvealtpyax(l,1) 

F4ACL VOL BAR (CL, Weel OCY (Pel )eaAqN 

FCT Mel se CT nbal eCY(SPat JOFEN( LT) MAIN 
ECLA UP IOCH  VeATE peal Ty NP) 

CANTINUE 

TFCNCALL EQ.tICaLE YS") 

ENOCHS) Foy 

CALL FSEMSCKrEVAEG, ALI SAIN, M21, d2™) 


CEFINE TRE MTAGSSSEN PRESSUPE 


DP ity Satyr 

ct 810 Lat,%s 

PCL oe JIaX(Tot,Jeony 
CINTINLE 


ald PRESSURE GRADTENT FRRACES TS vag ah vz 


e@ 120 Saiyts 

so" 120 Ts2,"1 

WAUCL eo JIBVRUCT Ye CPCI, JpeO(CILet esis (RmealR2(l)) 
28 938 Js2,'4 

“8 130 Sst,™3 

V2ZCTo SI BvZS(T SIO CPCI vc eP CL Jory ys (Rm*adle(d)) 
HE Tow 

Eno 


TIMEaIGe12812899 


09052000 
09033000 
0094000 
0999900 
0036000 
0097000 
0958000 
0039000 
oge0000 
0081000 
Coe2cod 
Oge3000 
9064000 
og8sooc 
oocedod 
0067000 
0088000 
0069000 
0070000 
60971000 
0072000 
0073000 
0074000 
9075000 
0076000 
90/7000 
0078000 
0079000 
0080000 
0081000 
0082000 
09063000 
0geudad 
00eSCcoG 
0086000 
oos?ogd 
008800G 
0969000 
oovouge 
eosTrcac 
covedge 
0993090 
0094000 
60950900 
00%8000 
0047000 
nesacgd 
oaVgcoc 
ayucdac 
OLeinge 
areeacyc 
C10 3cgo 
epnucoc 
O1uscyG 
0706000 
erercgd 


CakO IMOGE FILE EXZTAR(C FEN) om VESSE4N 05,29 Catewtusease2 TIMESS 4442912299 


wen “EMME OS™] 


te Staal ving Ah AQvICOO 
rzy PAGMar TE wBZ1,pag) Aguecad 
Ye DAMAMETER “LaMey  waameg,N{ Bhat,’ 2Bhe2 OoUusegd 
as CAGA*LTER “Pg %o1, APB EO} YQvelyl 
Se PARAMETEN MPL ae wef, PZ eMPed NPL aNPa] ,AP2ENPOD CovScue 
oe PARAMETER NAL ¥ az, ABLE Tah@Lkel 000g000 
Je REsle*® RCPRARINV RINVE PTTL A, ON MMS 0007000 
Re COMMA TAR ZICH) RSCMT) SORLCMT) pOPZ(M)SZICN) ,Z2CNL J eOZE CAL) ,OZACND OOVACOO 
Ga CUMIN SC VP SP LUV (MEQ HP2,NALK),RINVI (MPO, MPS,NALKS) PC OR (MO,3), vousagn 
100 t RTILt ACMP2),F (MP YNP) pNESIZECNALK) »TSCNALKS -SLMECNBLKD, 9020000 
ta 2 TE CNELED SE UL CMA) GE INCMPD) FAL CNPD FAM CNP) SAXCNP)ZAV(NP)s OOTIO00 
tee 3 SR CMB RP) CKCMP) SCY (NP) 0912060 
y3e CIMENS{(At. Cues vs (FP 2,2) 1913000 
{de ESUTVALENCE ESLER Y a ORING(N,S,SBLY)) Salacoo 
158 TE CY) SHASTZ2(1) e2 001§00C 
160 Be 9p NRe2,%ALK ANoletad 
170 TE (ME DSTE CMR f PONESIZE(*B) of 0017000 
14a Gq COATIMLE Coleogd 
198 re 96 wPat, ALK 9019000 
200 TS(™HOLISTE CAR Of aedéo0go 
ele Sg CANTINUE 00¢190¢ 
Zee pScrje@. 0042000 
e3n QP 1gd Trst,eP2 09423000 
ro of 140 Jat,3 agcadgd 
PL Ty e* 110 ta1,"P 0025000 
doa WCAR(T, 318000 0046000 
are tig CONTINLE goe7000 
rh RCSH(Itey,2.at,y 9028000 
eon NOgslecijet 0029000 
308 OS (30 Ji82,%9S 9030000 
sie C7 139 T82,"P3 0031990 
Sea ROORCE SPB Coae(pyeACIA (Tat, 2 ear (JL aFCOR (Tel eAB (Le Jie oo320q0 
33a ARCAMCT 2S veCX(T)aRCCR(Lot,2) ICV (1) 2053900 
34a L3E CANTISLE 00sucg0 
354 C7 140 Yat, “Pe 0055000 
Joe ACOR(T, 1) MRCAP (1,2) 0ose000 
37s ACOK(L,2)aaCsR(7,3) 0037000 
ye 14q CONTINUE 0058000 
39a 130 CIN IMLE 0059000 
00" of tus [et,“e2 0040000 
aie PINVCL Lg fet eRCOR( Lote) oo4taae 
42e OLeMVICTH,TIBRCPR(Loted) 6042000 
a3e 14u§ CANTINUE 0043000 
aaa 19§ CONTINUE 00440000 
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CONTINUE 
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of te? Jat,“© 
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CONTINUE 
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CONTINUE 

CANT INUE 
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END 


TIMEstust291E299 


vourcod 
Ocracos 
couseoo 
Coratot 
v04§eao 
90%eC00 
0097000 
0008000 
00%9000 
0010000 
gartega 
0g22000 
oo1scce 
oolacgd 
o0}$900 
OGtecod 
00717906 
votsevyd 
0019000 
0020000 
00€1000 
0022000 
9ad3000 
00¢4000 
90estg0 
00¢0000 
00¢70900 
00¢8000 
0029000 
0050000 
0031000 
0032000 
0053000 
09034090 
0035000 
0046000 
0037000 
0058000 
00359000 
0040090 
0041090 
oo4etod 
0043000 
o04ucge 
0045000 
0040000 
9047000 
0048000 
90499090 
0050000 
0031000 


CARO IMAGE FILE ECT TAR (CIEE) we VEPSTORN US,2% CATESINs 26/82 TImEm1 912812999 


ene “EMBE* HE*2 


te SLBRS TINE SGMOCK FRMAE SAN, AIN eo d21, heh) 0001000 
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oe FEAL OY ROMA RINVSRINVIZSTILOA 9006000 
Te CAMMAN EVP SAT UV (ME 26 MP Sg NALK) GQINV I (MPQ,MP2,NBLK1),ACOR(MP,3), 0007000 
i ae 1 RTTLL (MP2), F (MP, NP) WAST Z2(NBLK I, ISCNPLEK),SUXFCNBLED, 6008000 
j Ge 2 TE CRELE DoF TT CMP) SF INCMP) FST (NP) FAM(NP) AX (HPV ,AYCNP), 0009000 
100 3 AAims AP) o CHIMP) CY(NP) 0030000 
ite CIMENSTON x (HP, R8) 0011000 
' tae OF 90 NPw1, VALE 0012000 
| 130 Sunk (82890 0013004 
1 {ae 9q CANTINUE 0014090 
i 150 Ce 9S NPat,holh 0015000 
i tee en 9S [ao,nPt 0016000 
| 17s SUMF (NB) SSUMEC(NE) OABS(ECIZIE (' )01)) 0017000 
| 188 95 CANTINLE 0018000 
i 19¢ C7 Seq APset, ate 0019000 
20e TE (SU*FE (A) GT tt) Gt TS 96 Codo00d 
i ais Sieh (Neate o0¢1000 
4 ees Me CURTINUE 00e2000 
| ase NSTaRTal 0023000 
{ ge o2 199 TTYat,S 0024000 
2Se 99 206 NAwNSTART,AHL* 0025000 
t abe ISP¥alLS(nA) oy 0020000 
; 27a TE™2alE£ (NH) OQ 00e7T100 
{ 240 oe aus Jase, tere 0028000 
H 29a oO” 205 1202,%P1 0029000 
} 308 RCT e SOL BCE (Teen (Ti ax(Tel,JeaV(J)Ox(T Jel eGAlIeJ)e 0030000 
I Ste WC LoS IOCKELI aK (Tet edd IsCV CS? 0031000 
i 32a 206 Chr Tt ut ogs2ecgc 
B30 TE CNM SES MHLA) CR OTF 200 0933000 
jue o° See Itai,t) 005400 
; 350 JISTE CMR) of 0045000 
Joe ee 215 Ie2,'et 00560000 

370 ATLA CLTeL AKC E, SISLde CECT Si loA (Tax (Lotg dt ear(sije 0057000 | 
: yGn UeCLoeSlet pos (L si yexit, Ji) eCXCTIOXCT OL, JL /C%C 1) 9056000 
; 392 ats CTNFINLE 00 s900C 
i 408 Azavd,@ voan000 
i ate 1 2y¢ Tat, e2 1041000 
i uae asedQenaRScPTTL¢ acr)) 0042000 
fl USe ate CP TINUE rg43sco0 
| ade BSR AP/SLME (OA) 0046000 
' aSa tPCAS LEN tL) Ge 19 2k0 0045000 
i uoa 2" 297 Jai,3 0046N00 
: ule -* ay) Tei,<F 0ga7Uuod 
i" are erento l,Jjysn,n 6998000 
ute tr CONTIN LE 0949008 
d 506 7A 2d Jar," P2 0090000 
t gia PCN SON 2) 8U0y 9051000 
{ S2e oF 225 Jtar,”2 0952000 
y S30 ACHH(TOL,2BRCHE CJeLeZ)*QTILUACIEIMPIAVO JI, J,nNB) onszone 
| $40 = 223 CINTINUE 0034000 
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See O° 225 Ja2,"F1 00960000 
57s RCaM(J,1)80.0 090357000 
§%e G* 225 Ka2,“Pt 0998000 
§%e ACORN (S, 1) BREOOP(y,1) ORCOR CH, Z)ORINVI (Kel, Jet ghReoy) 0059000 
00s 228 CONTINUE 99000000 
ole C* 220 182,"P1 0061000 
eae KCTOTS(NB) CL, LS (rd) oRCOR(T, 1) 0002000 
Se = 226 CANTINLE 0063000 
ode 251 C*NTINUE Ngeacgd 
eSe CALL SS*ECKOTSCV AD, TECNA)) 0985000 
oe $22 CANT ISLE 00909000 
6?e = 230 C8MTINUE oge7C00 
oe J1aTE CNB) ot 0088000 
ede pe 220 Ia2,"et 0069000 
70a RCT Sires aCF CL Jt peak (LP ex( Lor ,JieaV (Ji ex(lpJlelje ooracad 
Tle WBA CL, Si eKCT St pec x Clyak(let,JivI7CV%Cdd) 0071900 
Tea 220 CONTINUE 0072000 
73a G01 CONTINUE 9a7soac 
740 200 CINTINUE 0g76000 
750 be 300 NB@tat, “apa 0o7S000 
Tee NPRNBLK eNBtol oo7eCgd 
778 TSP IBLSCNA) OY 0077000 
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201 CONTINUE 
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13e CARMMGN ZFOR SMELT, XTIME,TTIME, ISTEP, TSMC, ITAPE, TRY 0013000 
140 Cacti INnfune 0934000 
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16s REAUCS, LOadITIME 001000 
17s READ(S, 100) TER 0017000 
18e REAO(S, 100) TOUT 0038000 
1% REA0(S,109)3ISMO 0019000 
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ae oodsogo 
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280 Ge te 20 0026000 
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Jie ce 0051000 
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399 ¢ 0039000 
ade CALL UP Oo800g0 
ats¢ oo41000 
wae CAMBL Te ALL INVEISCIO TERMS 00420900 
asa ¢ 0093000 
aan Cayt 40VEC? 00448000 
a5e ¢ Oyesaud 
4oe ¢ COMPUTE VISCOUS TERMS 0988090 
ata ¢ 0o47000 
aa CALL OLFe 0046000 
age ¢ 0049000 
Soa ¢ MARCHING IN TIME 0030090 
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seec 0032000 
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5Se TFCUSTEP ,EG,1, AKO ITIVE EQ OVCEL Tag, e0EL? 0035000 
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gee c 0998000 
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ezs cc 0082000 
e3e CALL ChECK 0063000 
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oo ¢ 0980000 
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700 ¢ eor%oogs 
Tia €¢ WRITE KIQGTORY TAPE EVERY ITAPE STEPS 0071000 
Tae c oe72000 
730 TECMOOCI STEP TTAPE) EG O)WALTEC2)ITIME OATAL,Oataa,? Oa730G0 
jae 90 CONTINUE 0g7a000 
71$e ster eo7so00 
Toe ENO oa7ecgd 
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CARO IMAGE FILE ECITERC(CIFER) oe VERSION 05,29 CATESIO/26/602 


ane MEMBER INIT 


760 


NAONANA 


aaaan onan 


aAaann 


aan 


1¢ 


20 
30 


a9 


100 
110 


120 


130 


135 


1 
2 


TIMEalanigs32936 


SUBROUTINE INIT 
PARAMETER MBat,nuas 
PARAMETER Migtay ,P2SMa2 NI BNoj NZBNOg 
COMMON JONESVRACMAALISVTICMANID,BICME NE}, VRA(MANLD, 
VTQCH NED BAChE NID VRECM NL) DS VISCM NE), 
BIC NADP CMI NED, VZCMI AN) 
COMMON TWO 7A4 (HM), RACH1) SORE CML) »DROCM) ZICN) pZ2CNE) OZ (NL) OZ) 
COMMON STHA SRW, RHORCNY) sBV2(N) ALPHA, BENDA, BNDG CORY pGyhu (mM), ZK(N) 
COMMONSFORSOCELT,XTIME,TTIME, ISTEP, EMG, ITAPE, TOV 


INITIACIZE ALL CEPENCENT VARIABLES AND CONQTANTS 


ALPhA 39 THE NONOIMENSTONAL SMOOTHING COEF, 
FOR TIME SMHAOTHING IN SUBROUTINE FRRRO 


DEL Ta9Q0, 

aALPpMaso,10 

Ga980, 

Lare3o 

CORI 28? 2722EeSeSINCLAT a3, 14159/100,) 


DEFINE RADY AT GRID PAINTS AND ALL GRIC INTERVALS 


Of 10 1e1,™1 
0R1¢32820,€5 

R1(3) 80, 

00 20 182, 
RECTIBRAL (Tel Oh 1CI01) 
oe 30 81,41 

RACK BO.SeCRI CT ORICON?) 
ORQC1 wee (P24 eh1(1)) 
ORQ(MIa2 a(R (myeRacMi)) 
08 40 Iee,ni 
DRACTISBAI™CT IRQ Te1) 
MANSMAXMAG(OR1) 
ORMAXSOP1 (MAX) 


OEFINE ALL 02'S 


of 100 Jat,‘ 

Z1(J)8(Jot)*200,€2 

of 110 Jat,Nni 

OZ1CJIsZ1CJetdezi1¢J) 

OFAC I m2, 0(0,SacZzicryezt(1)) 
22(1)890,50022¢1) 

06 t20 Jaa,Nt 
OZ20JI80, Sel(O21( J) s0Z1¢Jet)) 

2265) 8Z22 (Jel) *O72(0) 

OZ2CNISE A CZICNJOZECNO1 00, SaOZICNE)) 


A ARO 6 ARE CONSTANT USED IN SUBROUTINE BeLNOV 
FaR CONSTANT OV, ANG VORT, CONOITIENS 


BNOABAL CHL) /RI CP) 
BNOGRCRICME DORI CM) MORE CME DACCRECME ORE (M2) ORI CM RORE (w2)) 
anobed, 


OEFINE CENSITY RELATED CONSTANTS 


RHOBL. 

06 130 Jai,N 
BVaCJI a! .£eb 

T3VB0. 

00 135 Jat, 
TOVSAMAXi(TAV,avacJ)) 
TOvVE1Le/SQRTCTAV) 


DEFINE PARTZENTAL AND VERTICAL DIFFLEION COEFFIENTS 


CPEPMEO COZAORE Ct et2/CELT? 

COEF ZEO,QO1FOZI CL) e*2/0ELT 

Of ta0 Tat," 

WK CT IOCOEPHE(S oF, aEXP(OFL OAT (Mie! )/7.)) 

06 1$0 Jai,Ni 

BKC IIOCHEPZE CL oS ACEXPCOFL OAT (J01)/5,) SEXP COFL GAT (NIi0J)/9,))) 
RETURN 

ENO 
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0001000 
cooeced 
0003000 
0ocecgd 
6005000 
0006090 
0007000 
0008000 
00veced 
0010000 
0011000 
0012000 
0013000 
001G000 
0015000 
0016000 
0017000 
0018000 
0029000 
0040000 
0021090 
00€2000 
00¢S000 
0024000 
0025000 
Vgeecge 
0047000 
00¢6000 
0029090 
0050000 
0031000 
oo32cq0 
0033000 
eoSeogo 
003Scge 
0050000 
0037000 
oodacgo 
0959000 
o4ocge 
oassoge 
0042000 
oo4scoo 
eosucgo 
0048000 
0090000 
ood 7000 
0048000 
0049000 
0090000 
0091000 
0052000 
0033000 
00Suoge 
009$000 
0099000 
0097000 
0056090 
0039000 
0060000 
0001000 
0082000 
0963000 
0gea000 
0085090 
0ob0000 
0087090 
0008000 
00849000 
0070000 
0073000 
0972000 
0073000 
0076000 
6079000 
0070000 
0077000 
coreced 


cet mein St 


CARO IMAGE FILE EDTTOR(CEFE*) eo VERGTEN 09,29 DATES10/26/82 TIMESLGsigss2ese 


aan MENGER START 


wy SUQROUTINE START 0001000 
ae PARAMETER “821,086! oov2cge 
3a PARAMETER Miata ,r 28%e2,NiaNot Nganeogd 0003000 
an COMMON SONESYRA CH ALIS VTICM ND DS ALCML NE) VRACM ANID, 0004000 
Sa 1 VTACH OND SASCML ONL oVRECMONE DG VISCMANS Dy 0005000 
oe 2 BICML NLP CMINET SVZCMI AN) 0000090 
Te COMMON s Tw 720 Cm) gH AZCH1) oDRECMI) SORZCM),ZICNIZZZCNE DD COZICNI),OZ2(N) COUTEGA 
Be COMMON sTHR INO RHER(NT) pBVZCN) ALPHA, ENDASBNDB CORT sGyhk(™),2K(N) 0008000 
Ge PARAMETER NOaZaranjoniany o00¢0e0 

100 OIMENSTON OATAL(NC) sOATAZ(NC) 0010000 

ite EQUIVALENCE (DaTAL,VRt)s (OATA2, VR2) 0011000 

fae c¢ 0012000 

1ye ¢ INITIALIZE “ASS PIELOS FOR A THEORETICAL QING 0013000 

14e ¢ 0014000 

150 riiat 0015000 

16a r1ae13 Ootecoa 

i?a T2ystl2er 0017000 
188 Tagan 0038000 
198 ar aGeU 0002 0019000 

200 of Lo falttstt2 0020000 

ate 16 BLCT ENE SAMAGMCOS (FLOAT CYeL 112/48, 23,12159) eG sRne 0021000 

22e Of 30 I8121,122 0042000 

23a Jo SECTONLIMOLCTIS,N1) SEP CoFLAAT (Lol2iol)/0,) 0023000 

aue Of 40 Jat,N2 0044000 

ase FACTSEXP (FLOAT (Jedi )/®,) 0045000 

aoe of Go Te1,9t 09026000 

a7a Oo BICl,J)eBicl 1) ee ace 0027090 

asec 9028000 

ase ¢ PRESSURE I§ GATAINED WYDROSTATICALLY FROP BUBYANCY 0029000 
Joe c 0030000 

3te c* So ren,st 0031000 

320 So PCI pi 2e0g, SeMhaed22(1 208141) 0032000 

330 08 69 Jag,ni 0033000 

340 of 60 Ist,m1 09034090 
3S» 60 PLZ oTIOPCT Jol og, SeRhOCOZZCJIA(BI CT SOI (I, Jol)? 0035000 j 

yee c 0g3ecgo 

Jlec TANGENTIAL VELSCITY I$ AY GRADIENT Balance 0037000 

Jee ¢ 0038000 

390 O08 70 Jat,nt 0039000 

aoe oe 70 fe2,mt 0040000 

aie PGFECP (IT, JIM (Tots JI /CUMAMORS(T)) 0041000 

ude RACBIO, SeCGRy any (J) )eaeQeni(l) aOGF 0942000 

aie JJaJ 0043000 

aun Itel 004u000 

use IF (HAU,LT.0.)60 TE 100 0045000 

doe Fo VV qCT ed meu eSeCORT ahi Cl )eSQPT (RAD) 09046000 

ava 0947000 

480 ¢ SET DATAQBOATAL FAR LEAPFROG 9048090 

avec 0049000 

50a ro 60 [wi,N0 0990000 

Sle Bo DATASCT adatal(y) 0091090 

S2e Cail eeLNoV 09032000 

S3s 9 90 Ya1,N0 0053000 

Gus 9q OATAICE Dedataacy) 0osd0c0 

§5e RETURN 0995090 

See 109 PRINT 119,21, JJ,PGF, RAD 0gSedcd 

S7e 11lo FORMAT(? RADICAL IN SLBRAUTINE START IG NEGATIVE AT (I,Jda',alSe 0037000 

$60 1' PGF, RAN Ot, 16 2612,5) 0958000 

LY) tev 0939000 


eae ENO voe0ce0 


CARD IMAGE FILE EOTTOR(CLFEN) wo VERSTON 05,29 OATESLO/20/82 TIMES1Gsigs32536 


eae MEMBER UP 


je SUBROUTINE UP 0001000 

ae PARAMETER MBQL,KB21 eocecce 

Be PARAMETER Hlamtoy, M2BMo2,NiaNeot ,N2ZBNO? 0003000 

ae COMMON SONE sVALCH ALIS VTECMONE) sBACMEGNI) pVRZCM ONL), TITTY 1 

Se t VISCH EAL pBSCMIANADAVRECHNLIZVTS(MANL), ooesooe 

6a 2 BCMA IP CME ONES VZOMI GND 0006600 

Ta COMMON /TMOsRE (Hy R2CHL) pORIECME) ORZCM) ZION) pZ2CN1) pOZL (HL) ,OZZCN) 0007000 

Be COMMON THR ARNO, RHERCNS) -BVZIN) p ALPHA, BNOA(GNOG CORT SG, HK CMI ,ZKCN) 0008000 

Ge PARAMETER MiNaw Tan 0009000 

100 DATA VZsF1N89,/ 0030000 

tia c 0032000 
ae ¢ BREQQURE 19 BRATAINED HYOROSTATICALLY FROM 8B 0032000 

430 ¢ 0033000 

140 08 10 Tat,mt 00348000 

15e La PCy pt) wed SeRHGe02261) e82(3e1) 0035040 x 4 
168 of 29 Jaa,nt 0016090 i 
17e 09 20 Ia3,M1 0037000 ‘ 
160 20 PCIe TIMP CT Jotpog, SeRmael7Z2(JIa(B2(I,3) e82(1,Je3)) 0038000 

198 ¢ 0019090 
z00e ¢ OLAGNSSE VERTICAL VELOCITY BY CONTINUITY eGUATIAN 0020000 
ate c¢ 0021090 
220 09 30 Jez,Nt 0022000 
ae oC 30 Ie1,m1 0023000 
24a JQ VZCTpXOVZCT, Jad oDZICIIACRACT OL aVAZCL ote TeLreR CT) aVRa(y,Jeol)) 0024000 
25a 1 ACORL(T 09, SacRiClor1) Mil) )) 0085090 
260 RETURN 0046000 

ara ENO 0027000 
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CARO IMAGE PILE EX TTAR(CIFEN) wo VERSTOAN 05,29 DATESIO/20/82 TIMEsIUsias32e36 


ean MEMBER BOUNDY 


{a 
ea 
Se 
aa 
Sa 
ee 
Ta 
Se 
90 
102 
ile 
120 
i3e 
148 
1Se 
Lae 
17e 


nNann 


SURROLTINE BAUNLY 
PARAMETER Magi, nal 
PARAMETER MiaMay  P2EMea,Nianat N2sned 
COMMON SONESVAL (CH ALI OVTLCMONE) SOL CHL ANID, VA2CM NID, 
1 VTSCM INGE BSCMLg NI) gVRSCMoNI DS VTISCM INI), 
BSCME NED cP CMEONL) »VZCMI ZN) 
CRMMONST HARE CM) ,RSCM1) DRI CML) SORZC4) ZION) pZZEN1) cDZ1(N1) -OZ2CN) 
COMPONSTHR/RHO,RMER(NT) sBVZ2(N) p ALPHA, ONDA SGNDB CORT Gy hi (HM), ZK(N) 


LATERAL @GUNDARY FOR TANGENTIAL AND PAOTAL VELSCITIES 
ASSUMING CONTINUCUS VORTICITY AND DS VERGENCE 


of to Jai,Ni 

VRa(Me JI SBNOARVRZ (MI pd) OBNOBR(RI (ML) AVRA(ME DS IOAL (M2) aVR2(M2—J)) 
10 VTaCMod @BNOAaVTA(ME eS) OONOBR(RICMSAVTSCML eo TIGRE CMA) AVIS (M2,5)) 

RETURN 

END 


oguticoe 
Oguecga 
0003000 
oooatco 
oouscgo 
ogvecegd 
0u07000 
9048000 
oeoeogo 
0010000 
0011006 
0032000 
0013000 
00184000 
0015000 
0036000 
0017000 


CAML IMAGE FILE 


eae MEMBER OFF 


eT TESTER eRe ——— 
- 
oOo eareWowhe 
@aenprnenpnennvesee 


$9e 
000 


o3e 
ode 


64a 


ECTTER(CIFER) om VERSIEN 05,29 DATER10/20/602 TIMEMLAgios32a3e 
QUBNSUTINE DIFF 0005000 ‘ 
c 0oczcgo \ 
¢ CamPLTE THE DIFFUSION TERNS 0003000 
c 090ec00 
PARAMETER M@Qt,neet 0005000 ; 
PARAMETER Mighes ,h28Me2,NiBNol ,N2aNe2 900ec0d : 
COMMON SONESVAL CoAT GS VTLCM ND) SBLCML ONE), VR2(MyNEDy oocTcad 
1 VTSCMONEDSBZCMI GND Dp VRECMONI) DVISIMYNE) > 0008000 
BSCML NL) SPCML ENED» VZCME LN) 0000090 


aan 


nnn 


COMMONS TWA SRE CM) RACHID) ORE CHL) SORZCM) ZL CNIS ZZCNS) SOZL CNG, OZZCND 0030000 
COMMONsSTHRARHG KHOR CNT) pOV2ON) ALPHA, BNOASBNOB CORT Goku (H),ZK(N} 0033000 


DIMENSION VROM Hid, VT CMeNS) pO CHIANG) 0012000 
EQUIVALENCE CYR VEL o CVT VTE), (8,88) 0033060 
0014000 

MORIZONTAL OLFFUSION GF RAOTAL VEL@CITY 0015900 : 

0036000 : 

DO 10 Jai,ni 0017000 : 
08 ty f#2,m1 0018000 
Lo VASCIoJISVR3CT JI oMK (TI MCC CVRCL OLS SIOVROT, JI) /ORICI) 0019000 
1 PCVRC Ze TOUR Clegg SI)JORI CIOL )SORSCT IVR CT eJISCMECEIERICI)) 9020000 
2 ODeSAC(VECL OL, Jp avRAC Te JII/CORICTMR2(T)) 0021000 
3 OCVACT SOWA (Teg yd) (CORE CT 01) oR2(103)))) og22cod 
0023000 
HORIZONTAL DIFFUSION GF TANGENTIAL VELOCITY 0024000 
0025000 
0% 20 J381,N1 0026000 
be 20 182,81 0027090 
ZO VISCI oe SISVTT(L JI OMK CTI OC ECVE CL ele JdOVT CT oJ I/ORI IS) 0028000 
1 SEVTC Tp TISVT (Lote dT) /ORUCIaqAORSCTIOVICT se JISCAS(TIMRICTI) 0029000 
2 OOeSCOVICT OL Spmv PCL, JI /COR1CIIMRACT)) 00950090 
3 OCVTCT JeVTCLer, J) ACORL (Tol) aR2(I01)))) 0031000 
oo3Z000 
HORIZONTAL OIFFUSION OF B 00330y0 
oosucoo 
08 60 J81,N! 0035000 
DA 60 1eg,m2 0036000 
00 GSCI gS BRRCI gt xome CVI a(CCB(Lol eS waCI oJ) /ORQCI 41) 00579090 
{ @CACT i ddeOCLete J) )/0R2(1) 2 /DR1 (1) 0038000 
2 FOeSOC(R(L Ole SIOACT oJ) I/CORSCT OL ARIA CT O$)) n939000 
3 HCBCT JIB Clore JIVACORECT eat CL3I33 0940000 
06 70 Ja1,Ni 0091000 
asi, seasct, JyehM (1) 0 (CBC SeII 0B 01,3227 (OR202)#OR1(1)) 0942000 
00058 (E (29d) 08 (1,4) 27 (DR2C2)#R1(2))) 9043000 
"68 By J8L,Ni 09084000 
ag Bee ee UME see 0045000 
mCBECME, JICB(MAsJIIZ(NRAC41 DARE CMI))) 0098000 
00470900 
VERTICAL OIFFUSION OF RACTAL VELSCITY 0048000 
0049000 
oe 90 Jsa,n2 00$0000 
of 90 [8a,m1 0051900 
be VAMC eo JIaVROCT,SVOZKCIIMCCVRCI Sot vovA (led) /022 (Jel) 60392000 
@CVACT  SIOVACT Sot) 702263) 702860) 0093000 
— 100 fa2,%1 0094000 
10G VAZCT ot IavASCL,FdeZKCUPMCVRC So 2IeVRCL41))/(00Z2062) 0021 (1)) 0035000 
06 110 Tea," 0090000 
WLG VASCTONDISVASCT NE DOZR CME DMCOVRCT NI DOVRITON2)/(CZZCNI)MCZLCNED) 0097090 
0098000 
VERTICAL DIFFUSION QF TANGENTIAL VELOCITY 0099000 
g0e0dg0 
09 120 Ja2,%2 0061000 
06 129 I92,"1 0062000 
120 VISCTeSISVTS(L J peZKCSIOCCVPCE Jat pour (led) 022 (Je!) 0063090 
1 ecvT (ly JeVTCLe Jot) ) 702265) 9/021 00) Ooeuogd 
08 130 te2,%1 0065090 
UFO VISCIDIOVTIT CTP OZRCVIOCVI CT» 2I0VIC3 4199702202) 00721 ¢61)) 0000000 
08 140 t82,%4 0967000 
199 VTSCEoNEOVISCENIIOZKCNE MCOVITCE NE DOV TCT eN22IA(COZ2CNE ) ONZICNI)) 00ba000 
0969060 
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CARO TmaGE FILE ECTTOR(CIFEX) es VERSION 05,29 DATFelds2es02 TIMEMIGstes32e30 


ane MEMBER OLFF 


700 
Vie 
yas 
738 
yas 
1$0 
768 
778 
780 
7908 
20a 
ale 
eae 


Onn 


VERTICAL OFFFUSTION OF 8 


oe 169 Js2,N2 
Of 160 Ilei,*1 
160 Pa JBASCT Jp oZM CSC CBICT Jol oaci, J) )/022( Jet) 
(BCI pd eB Le det) ) 7028099) 702163) 
ar 170 Yai,™t 
170 a3¢1,})003(I LpoZMCLIA(ACT, 2) 9863 xe) 
o7 180 1st,"t es Si suena eae 
a N 
149 Be CT ee LIMES CLONED ES ENCE GT CE TOR Cha Nets CRENNAPUE TE 
END 


0070000 
OoTLege 
0o72090 
Oo73C00 
0074000 
0075000 
0076000 
0077000 
0076000 
0079000 
0980000 
9061000 
0062000 


CARD IMAGE FILE ECITER(CIFEM) eo VERSTOR 08,29 OATER10/20/62 TIrM€miagiars2e3e 


ane MEMBER FRWRD 


ta SURRGQUTINE FRawat oguidgd 
ae PARAMETER wa2i,neg! 0002000 
Je PARAMETER Migtay,hZahe2,Niane] N2Bheg 0003000 
Qe PARAMETER NOaZavaht{omieany oousu000 
Sa COMMONSONE JOATAS (AD) CATAZ(NO) ,CATASCND) P(ME SNL) VZCMIN) 0005000 
os COMMON STHRARWO, RKER(NT) ¢BVZCN) ALPHA, BNOA,BNOB, CORT »G, M(H) ,ZK(N) 0006000 
Te COMMON SFORSOELT XTIME S TTIMES ISTEP, SMO, 2 TAPE, THY 0007090 
aec 0006090 
9a C REPLACE CATAS ITH THE NEw VALUES 0009000 
toe c 9010000 
1is 00 10 I8i,N0 0011000 
12a 10 OATASCEaOaTalr (yes, eCELTeQaragi(s) 0012000 
ise c¢ 0013000 
14a ¢ TIME SMCOTRING ogLacge 
1s¢¢ 0015000 
1oe IF(MOUCESTEP, 190), NE,0)G0 TO 30 0016000 
17 Of 20 I81,Nn 0037000 
188 20 OAVASCT MOATASCT IO COATAICT)OOATASCE a2, eOATAQ(T) OaLPwa 00180900 
198 Jo CONTINUE 0019000 
zon c 0020000 
aie c FRRWARO MARCHING 0021000 
a2ec ocd2000 
aia DO So Js3,N0 00230900 
ean Og OATALCT eNaTaacy) 00de000 
250 OC 50 Iat,Ne 0025000 
260 So O*TAlCT  ehaTas(y) 00¢9000 
a7 c 0027000 
28e ¢ ZEA ALT Catads FOR NEXT STEP 0026000 
age ¢ 0029000 
308 08 60 Iay,nc 0030000 
Jin 60 DATAZII) 80, 0031000 
320 RE TURK oos20c0 
330 ENO 0933000 
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+. rane vervepnnnnier 


CARU [maGe FILE ECLTER(CIFEM) em VERSTOR 05,29 D&TEatosaess2 TIreEstasiot3arse 


ann MEMBER CrRECK 


le 
ae 
Je 
te 
Sa 
be 
Ta 
Se 
qe 
100 
ile 
12e 
t3e 
jae 
1Sa 
16a 
178 
164 
190 
208 
aie 
22ee 
230 


20 


10 


190 


SLBNOLTINE CmECS 

DAQAMETER uO2!],nOa! 

DARAMETER Mistoy , »2EM%odsNiOhal ,\28Ne2 

CHHMEN SANE SVRL (CH ALISVTIC™ NS), BICME NL), VRA(MLNG), 
VTSCHALISASCMELNED, VRSCMNLI,VTS(MNL), 
BsCmy Ard sP CME NLD eVZ(ME AN) 

CHMMONs Te IsR 1 (H) BACHE) ORL CML) ORS), 21 (ND, Z2CNI) POZLCNI) pOZ2CN) 

COmM“ON FOR SOELT, KILME,ITIME, ISTES, LSMP, I TAPE, TaV 

OIMENSTON wtRK yh), MORK 2(N) 

oF te J81,N1 

ce 20 184,” 

MORAL CI OMRZ(T) sah OXi (1, ,VR2(1,5)) 

MINSMINMAG (ORAKY Ot 

CTenOwni (MIN) 09.9 

OTaAMINICOT/CELT) 

CONTINVE 

OTeAMINICOT, TAY) 

TF COT SGE PELTIRETLRN 

DEL T#O, Se0ELT 

PRINT 100,0ELT 

FBRMATC As /s/s,'eeoneeeDELT 1S CHANGED TO, IPEL1,2,!' Saneenteen!) 

RETURN 

EXC 
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0001090 
0002000 
0003000 
0008000 
0909000 
0990000 
0007090 
0008000 
0009000 
0010000 
0014000 
oot2000 
00135000 
0a14dgd 
00359000 
0016000 
00170900 
0018000 
0039000 
9020000 
0021000 
ogd2cod 
Codsco0 


ee nn Re casei eta tall ae eee A a Met a 


Area aa a 


CARO IMAGE FILE EL LITER (CIFEN) em VERGTAN 05,29 DATEBIOs20/82 TIMEDIO3$ 69392936 


een MEMBER AQVECT 


te SUBROLTINE AQVECT ocv1ogd 
ee 2902090 
Sac CerPLTE THE ADVECTIVE TERMS 0003000 
gee ocouggd 
Sa PARAMETER maat,ragl 000sdgd 
ee PARAMETER Miateal, »geMog,niatet ,N2aneg 00ug0u0 
te CHMMONSONESVRI (DED A VTICM ONG) ,BL (MINE), VR2(%,NI), 00070900 
fe 1 VTQCH NL), BS(% Ly W1) p VAS (MON) pVTSEM NL) y 0008000 
Ge 2 BSC*r Mad eP CME NED, VZCMIEAN) 0909000 
108 COMMON sTRASZRE (CH) RACHEL) GORE CME) ORIACM) ZICN) ¢ZZ20IN1 9,029 ¢NI),0Z2¢0N) 0010000 
ile COMMONS THRsAWA RRERINE) pAVZ(N) ALOHA, BNDA,BNDB CORT, GmMuCM),ZE(N) 0011000 
lee DIMENSTON VROM,AT) VI CMANT) SOCMELNY) ooi2cga 
13a EQUIVALENCE (VR, VG2Ie (VT, VTS), (8,82) 00133000 
‘ 14a ¢ 0¢14000 
i 199 ¢ M@BIZONTAL ADVECTION FOR RADIAL VELOCITY 0015000 
i tee € Oortecgd 
ifs 06 16 Jat,st 0017000 
180 08 10 Isa,et 0018000 
190 10 VRECT eS Bed 28a (VAC Le JI OVACTote J) eCVACS pJICVACTore J) CRI CTel) ootscasd 
200 t OCVACE OLS TI*VACT  JIIMCVACT Ole SIOVRED IIIZORICI5) 6020000 
ate 2 #vR3CI,J) 00421000 
a2e¢ 00422000 
ayac WORIZONTAL ADVECTION POR TANGENTIAL VELOCITY 00@30G0 
@4e¢ 002a000 
ase 08 20 Jai,Nn! 00eS000 
ace 08 20 Le2,m"1 0odeceo 
ave 20 Te uae dds td pale UL ie ties be SyevT Cor, JI) sORL (ket? 9027000 
ase CVA CT ets JIOVRIT JIINCVI (TOL, JIOVI CT oJ? /0RI(3)) 0026000 
298 : ovT3(3,J) 0029000 
joa MORTZONTAL ADVECTIAN BER BUCYANCY 09050000 
yee 0031000 
Jae 08 60 Jat,n! 0032000 
330 Oe oy Ilaa,m2 0033000 
348 60 seeth JBOSS Cis TIO SHCVACI A JIC(BCL oe JIeh (Tete) }/0R2(5) 0030200 
SSe VRCLot pJ)@CB (Lote J) @BlT, 2) )/0R2164)) 0035000 
Joe "56 70 Jei,ni oodecce 
370 70 B31, IIMAS C1 J) 00 ySAVA (2,3) 8(B62y S986 Led) )/0RZC2) 0037000 { 
t jae co 80 Ja1,Nt 0038000 
‘ 398 BO BIC*L eT IABIEMI VO SAVRIMG pJIOCBCM1, cI OB(M2,J)) /OR2(HL) 0939000 
7 ace c 0040000 
alec VERTICAL AOVECTION FOR RADTAL VELSCITY 0041000 
ade cg 0042000 
a3e £9 90 Jee,Nn2 0043000 
ade 06 90 Ima,mi 0044000 
aSa 90 VRZCT oe SIMVAS(T Jade dSaC (VEC Lede SIOVZC Lo TIIMCVAC Te J IGVR OE Jet)? 0945000 
aee { POZACIIGOCVZCT Jos vovZ (Jos, Jos IacvR (Ip Joi ovRlyes)) 0046000 
aye 2 70220 501)) 0047000 
abe 05 9S Jea,mt 0028090 
ade 9S VASCT pL IBVAS(T 1 pede dSal(VZC Pe 2)ovZ(lot, 2) ACVRCT,2.e°VR( 7 51)) 00499000 
g0e { AOZ%2C2) 0030000 
Sis 08 % Tea,mt 0091000 
S2e Ge VRICT ONS SVASCE NE OO asa CVZCLei Ni dovZC To NI)IMCVACI NI JeVRCT,N2)) 8092000 
$30 1 70Z20N13 0993000 
$4e¢ 0954000 
: B50 ¢ VERTICAL ACVECTION FOR TANGENTIAL VELSCITY 0935000 
: See C 0036000 
: $76 08 890 Je2,N2 0037000 
i $80 p6 :00 sa2,"1 0998000 
H 990 10Q VISCL eT VTS CE JO kSaCCVACLop,JIoVZC Lod ACVT CL eS IoVI (I, Jol)) 0oSecqd 
i eda 1 JOZSCIIONVZCT Jot yevZ (Let, Jol Pecvie(IedJeolyeovil (rsd)? 09¢0000 
ole 2 AOZACJes)) 00619900 
e2s Of 105 Ie2,"t 0062000 
o3e 108 Gee Peet RANT) Cah vems da eRe Cig EIS VERT ENR DOE Sete chenae 00963000 
ode 7022(¢2) 0900000 
e3s ‘0 toe Te2,"1 0905000 
ebe 106 VISCTomtvevtaer, NEDO SSACVACI ND IOVZCLOLINEDIOCVICI ANE OVI (I gN2)) Gaoedad 
eve 1 702204) 0007000 
e8e ¢ 00b8090 
eo ¢ VERTICAL ADVECTION FOR 8 0069000 
700 ¢ 0070000 
Tie 68 190 Je2,%2 0071000 
Jas 00 160 1e1,"14 0072000 
130 180 per SMBS CT, JpeOSACVZCT JIMCBCI Jeb (I ,Joi})/0Z2 (J) @o73000 
Ie ove (ty Jot pacLpdote8 (yd) )/022901)) 0078000 
730 ‘oe 1e0 Tei,"3 007S$000 
tee §=6160 BCL» V)MOS(T, 00 ,SAVZCL,2)e(M|(T,2veBatty1)I/02ala) 0678000 
778 OC 170 Tet,"1 0077000 
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CARO IMAGE FILE ECTTAR(CIFER) wo VERSION 05,29 DATERtO/20/22 TIMEslarlas32236 


: ate MEMBER ADVECT 
| 788 170 BICTONUDGBSCI Ms dee SOVZCIT Nt aCOCIoN1 ped (TeN2) I ZCZ2CN1y vo7eCo0 
798 C 0079000 
1 ede ¢ INERTIA TERMS FOR HORIZONTAL MOMENTYP 0980000 
i sie 0081090 
ade OF 190 Jetynt 0062000 
63a of 110 {m2,"1 oo83Coe 
60 VAZCIo TI OVASCITIOVICT JI MCVTCE eT) /RA CT OCORY) 0088090 
! ase ULG VISCTeSIMVTSCL JIOVACT JIOCVICTs SI SPICITIOCORT) 0085000 
; eee € 0080000 
are PRESSURE GRAOTENT FORCE 0087000 
e8e € 00880y0 
ade OO 120 Jat,%s 0969000 
qQ0e oe 120 Is2,"1 00900900 
10 £20 VASCIe JIBVAS CL JdeCPCI, Jp oP (Tete JI) /(RHOROR2Z(1)) 0091000 
G2a¢ 0092090 
O30 C¢ STRATIFICATION TERM 0043coq 
Qua ce 0044090 
950 08 130 Jat,%4 0095000 
qd6e os 130 Tat," 09%e0y0 
g7a USO GSCLASBAS(T, Jeo SeCvICT, J aBvaCJev2(leJolpeavacjei)) 0097000 
Qae RETURN 0096600 
99 ENC 0099000 


5 ns tang An snail n> 
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nr oS Preteen teppei 


CARO [MAGE FILE ECT TGR(CIFEX) wo VERSION 05,29 DATERI9/20/02 TIMESLasses3as30 


een MEMBER PUTOUT 


te SUBROUTINE PUTALT 90014006 
aa PARAMETER MO2L, NEE! dovu2gdqd 
3a PARAMETER Hi shat, h2EMe2,Nighol ,Naened 0003090 
Ge COMMON JONESVAL CH ONT DG VTECH ANID ZALCMEGA EDA VRELM ONT), oocudgs 
Se { VTQCM ND eBSCMD NTI oVASCMEND) PVTSCMONIL Dy 0005000 
be 2 BSEML,NTDeP CMS NED,VZCMI,N) 0000600 
ve COMMON STRASAL (My pR2CM1) PORE CME) PURZCM) ZION Gp Z2CNI) POZA CHI OZACN) OCUTCOS 
Ae COMMONS THR ZRHE RMERONT) pAVZCN) pALPWAENOMSBNOD CORT sGyhe(™),ZK(N) 09008000 
Ge COMMON SFARSOELT KTIME,LTIMES ISTEP, SMO, J TAPE, TBV 0009000 
19e ¢ : 0010000 
ye ThIg SUBRGUTINE PRINY BLT FIELOS FOR A QurtK Lean oolicgo 
see 0012000 
13e OIMENSZON [CUM(m,n) 0013000 
tae 700 FORMATO As//¢" wACTAL VELSCITY (CM/S) aT Tet, lo,* mw) 0014a000 
15* 70§ FOAMAT(As/74° TANGENTIAL VELOCITY (CMsS) AT Te',ley! Ht) 0038000 
poe TiO FCRMAT Ce ssy* VERTICOL VELECITY (CM/S) AT Tet, Tey! HT) e0teooo 
176 T16 FORMAT(/sss—¢* BUOYANCY FIELO (8,061) AT Tat, le,* &') 0037000 
168 725 FORMAT(/s77q' PRESSLBE (910 DYNEsCMeeg) AT TIMER, 16,* &*) 0oteona 
198 720 FARMATCIMI,/s//,20Ks  orerorere OLTPUT AT TIME Bt,te,t ORI, 0019000 
20e 1 FRede' CY ISTEP at tte? Farvereceoesool) ooz20000 
ais OSVOxXT IME s 86400, %0,0061 0021000 
aga ORINT T20,TTIME,CAYe LSTEP 09¢2000 
23s 06 tu J81_NI 00¢3000 
aan oe to ret,M 0024000 
25a 10 LOUN(I,J)evRact,J) 0025000 
eee PRINT 700,I TIME 00¢0000 
eve CALL MAP CIDUM,R1,227%)N3) 0047000 
260 Of 20 Jet,N} 0028000 
aga Of 20 181, 0029000 
308 20 1OUM(lL J evTacted) 0030000 
ste PRINT 7O0SeITIME 0051000 
32s CALL MAPCIDUM,RE,22,M,N1) cos2uoe 
330 09 30 J81,N 00330c0 
jae co 3o ret,mt 003000 
35* $o 10UN(L,Jpavic(t,J) 0035000 
3oe PRINT 71O,1TIME 0038090 
37s CALL PAPCTOUM,RE,Z1eM1,%) oo37ou0 
yee Of Yo Jat,Nt oadacgo 
398 On 99 [a1,Mt 0039000 
aon 4g 1CUM(1,I)#a2(T 9981 E3 0049090 
aia PRINT TIS,ITIME 0041090 
wea CALL MAP(TOUM SRA, Z2e%1,%1) 0042000 
aga c7 SO Jai,Nnt 043090 
aan Of So Jatyet 0osaugd 
uSa SO yDuUM(T,J) ae ly, J a1,Es1 0045000 
aoe PRINT 725, 1TIME 0C4eg0ad 
uqe CALL MAPCIDUM Reg Z2eMt,%1) 0047000 
aha RETURN 00460900 


aoe END 60490900 


CARO IMAGE FILE ECSTORICTFE™) om VEPSTOR 05,29 CATERINs20783 TIPEataster 32136 


eee MEMGEN “AP 


1s SUBRSUTINE “ah lage lg ¥ hy NA) 9eulegd 
ee PARAMETER “m2i,neel 0002090 
3a DIMENSION ACM), ZENA) 0ovso00 
ae INTEGER A (MAND, TE (MDG TZ(N) oguscge 
Sa To FOR™AT CIES, 7x,2818) 000S000 
be AQ FORMAT (IMS, Tu, 3x, 2518) 000e0u0 
Ta MPaM[NQC 85, ¥4) vovroyc 
Ga AA LO Ta1,MP 0908000 
Ga fo IRC IMRCT) al £esec,! vgveugd 
108 ve 20 Jei,NN egto000 
tte Zo T2(J a2 I) at Eeceoc,! 0011000 
12e PRINT 70 9912000 
tie PRINT OCIA CII, T81,"F) vargsogo 
Lue PRINT 70 0914000 
1Se 08 30 JJai,Nn 0015000 
lee JanNeleds oo1e000 
17e 3o PRINT BO, TZCJI, (ACT eS), Tet YP) 00170¢0 
188 RETURN 0016000 
190 END 0019000 
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